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Abstract 

We compute the nuclear force in a holographic model of QCD on the basis of a 
D4-D8 brane configuration in type IIA string theory. The repulsive core of nucleons is 
important in nuclear physics, but its origin has not been well understood in strongly 
coupled QCD. We find that the string theory via gauge/string duality deduces this 
repulsive core at a short distance between nucleons. Since baryons in the model are 
realized as solitons given by Yang-Mills instanton configuration on flavor D8-branes, 
ADHM construction of two instantons probes well the nucleon interaction at short 
scale, which provides the nuclear force quantitatively. We obtain a central force, as 
well as a tensor force, which is strongly repulsive as suggested in experiments and 
lattice results. In particular, the nucleon-nucleon potential V{r) (as a function of the 
distance) scales as r~^, which is peculiar to the holographic model. We compare our 
results with the one-boson exchange model using the nucleon-nucleon-meson coupling 
obtained in our previous paper.^l' 
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§1. Introduction 

Nuclear force, the force between nucleons, exhibits a repulsive core of nucleons at short 
distances. This repulsive core is quite important for large varieties of physics of nuclei 
and nuclear matter. For example, the well-known presence of nuclear saturation density 
is essentially due to this repulsive core. However, from the viewpoint of strongly coupled 
QCD, the physical origin of this repulsive core has not been well understood. Despite the long 
history of the problem, it was rather recentf^' that lattice QCD could reach the problem]^ 
and of course, any understanding of it based on analytic computations is quite helpful for 
revealing the basic nature of nuclear and hadron physics. 



The recent rapid progress in applying gauge/string dualitj^^^ to QCD, holographic 
QCD, has been surprising. Now, it has been made possible to compute various observables 
in hadron physics such as spectra of mesons/baryons/glueballs and the interactions among 
them. Although most of the works rely on the supergravity approximation that works for 
large Nc and large 't Hooft coupling A, it turned out that the holographic QCD reproduces 
quite well the properties of hadrons not only qualitatively but also quantitatively. 

We app ly this gauge/string duality to the problem of nuclear force. In our previous 
we computed nucleon-nucleon-meson couplings, using the holographic QCD on 



paper's! 



the basis of a D4-D8 brane configuration in type IIA string theory,'^''^^ which incorporates 
chiral quark dynamics. This amounts in principle to computing the large distance behavior 
of nuclear force, given that the potential between two nucleons can be understood as an 
exchange of mesons between them. In this paper, we take one step further: by directly 
solving the two-nucleon system in the D4-D8 model of the holographic QCD, we find a short 
distance scale of the nuclear force. In fact, we find the repulsive core of nucleons. 

First, let us briefly summarize what has been computed for baryons in the D4-D8 model 
of the holographic QCD. The D4-D8 modeP''^ of the holographic QCD describes a strong 
coupling regime of massless QCD at low energy, in the large Nc limit with large 't Hooft 
coupling A, for a fixed number Nj of flavorsO The low-energy degrees of freedom on the 
flavor D8-branes in the holographic geometry of Ref . |T2j) , which are basically the Yang-Mills 
(YM) fields in five-dimensional curved space-time, give Kaluza-Klein towers of mesons, while 
instantons in the YM theory correspond to baryons in low-energy QCD® (this is based on 
the baryon vertices in gauge/string duality'^''^ and the fact that branes inside branes are 



*•* See also Ref. |2| for a study of the interactfons between nucleons and hyperons using lattice QCD. 
For a review on the theoretical aspects of nuclear force including that of short distance, see for example 
Ref. 13). 

See also Refs. [T11),[1T])-[251) for closely related works. 
Introducing massive quarks in the model has been discussed in Ref. [TT|) . 
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represented by solitonic instantonJ^). Here, in our terminology, the instanton is a gauge 
configuration that is locahzed in spatial four dimensions in the five-dimensional space-time. 
The baryon number is identified as the instanton number in four-dimensional space. Since it 
is conserved in the time direction and localized in the spatial directions, it behaves as point 
particles that are interpreted as baryons. Quantization of a single instanton a la moduli 
space approximatiorJ^''^ gives rise to a spectrum of baryons including nucleons.'^ In our 
previous paper,'^'' we computed the static properties of the baryons by evaluating the chiral 
currents in the presence of the instanton: charge radius, magnetic moments, form factors 
etc. were computed, in addition to the nucleon-nucleon-meson couplingso This analysis 
is reminiscent of that by Adkins et. al.'^ for Skyrmions.^**^! In fact, the holographic 
descriptio n ra imics the relation between the Skyrmion and instantons found by Atiyah and 
Mantoij ^***^! . The physics of finite baryon density and nuclear matter has been explored in 



many papers recently, and we do not describe them in detail here. 

Next, we briefly outline our method. The one-instanton analysis given in Refs. [TSjl . fT^ 
revealed that the desired configuration with instanton number 1 can be obtained simply 
by considering corrections to the BPST instanton in four- dimensional flat space.'^ The 
corrections are due to (i) overall U{1) part of the YM gauge fields coupled to the instanton 
density, and (ii) curved space-time along the extra dimension in the five- dimensional 
space-time. These corrections induce a small potential in the instanton moduli space, fix 
the size of the instanton to be of order 1/(a/AMkk) (where Mkk is the only parameter 
with mass dimension and gives the meson mass scale), and give the quantization of the 
instanton in the moduli space approximation. This type of analysis can be extended to 
our case of two baryons. If the two baryons sit close to each other so that the distance r 
satisfies r < (9(1/Mkk), we can use two-instanton configuration in the fiat space as a starting 
point, since the effects of the curved space are small. The properties of the two-instanton 
moduli space are known, concerning not only its construction via renowned ADHM (Atiyah- 
Drinfeld-Hitchin-Manin) method,'^-' but also the metric in its moduli space (see Refs. [32|1 - [36|1 
for some of the papers relevant to our computations). We use them explicitly as a basis in 
a manner similar to the one-instanton case, to explore the physics of the nuclear force, i.e., 
the interaction between two baryons sitting close to each other. 

We compute the additional potential induced in the moduli space, due to the presence 



By using five-dimensional spinor fields introduced as nucleon fields on the D8-brane, in Refs. [T ^ .[20 1 , 
the static quantities of baryons were computed. See also Refs. I ^ .l ^ . 

The analysis of the Skyrmions based on the four-dimensional meson effective action derived from the 
D4-D8 model is given in Ref. [?7|. 

***) Xo describe the nuclear force, a two-instanton configuration was used for this Atiyah-Manton ansatz 
for Skyrmions (see for example Ref. 
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of the two instantons. The analytic form of the moduh Lagrangian can be obtained in the 
asymptotic expansion of r. This includes corrections to the kinetic term, coming from the 
metric of two-instanton moduli space. Specifying the two-nucleon states by tensor product of 
single-baryon states obtained in Ref . flSjl . we can evaluate the nuclear force for given nucleon 
states. 

Note that this "asymptotics" means a large distance in the region r < 0{1/Mj<i^). There- 
fore, in the standard terminology for the nuclear force, our result is for short distances. In 
addition, we use the asymptotic expansion in r, so our analytic formula of the nucleon- 
nucleon potential is not for nucleons on top of each other. However, this is sufficient for 
seeing the repulsive core of the nucleons. 

We find that our final expression for the nucleon-nucleon potential, f l4-45p . is repulsive, 
and has dependence. This r-dependence is peculiar to the four- dimensional space, not 
the three-dimensional harmonic potential. The appearance of the potential is due to the 
extra holographic dimension, thus typical in holographic description. Physically speaking, 
the Kaluza-Klein summation of all the meson states in the tower produces this new behavior. 

The main reason why the force is repulsive is that the instantons carry electric charge 
of the overall U{1) part of the YM fields on the D8-branes. This electric charge is supplied 
by a Chern-Simons (CS) coupling on the Nf D8-branes, and is nothing but the baryon 
number. The U{1) force is repulsive since the instantons have the U{1) charge of the same 
sign. There are some contributions from the SU (2) gauge field that give attractive potential, 
but this SU{2) force turns out not to be strong enough to cancel the U{1) repulsive force. 
The Kaluza-Klein decomposition of the U{1) part of the gauge fields provides a mass tower 
starting with u meson as the lightest vector meson,'^^ and so, our computation shows that 
the repulsive force is partly due to the u meson exchange. Not only the u meson but also the 
whole massive mesons participate in the nuclear force, and as a result, the nucleon-nucleon 
potential becomes 

One might wonder whether it is reasonable to sum up the contributions from all the 
massive mesons, since the model deviates from QCD at the energy scale higher than Mkk- 
However, there are some lines of evidence suggesting that the results obtained by summing 
up the infinite tower of massive mesons are better than those obtained by only taking into 
account the first few modes. For example, in our previous paper, we showed that the 
electromagnetic form factors for the nucleon are very close to the dipole profile observed 
in the experiment. This result is obtained by summing up the contributions from all the 
massive vector mesons. If we only take into account the rho meson, the form factors can 
never be close to the dipole profile. 

The organization of this paper is as follows. First, in ^ we describe our strategy. 
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together with a brief review of the instantons in the D4-D8 modeL In ^ we obtain an 
effective Hamiltonian for moduh parameters of the two instantons in the modeL In §U using 
the wave functions for nucleon states, we evaluate the nucleon nucleon interaction potential. 
We decompose it to a central force and a tensor force. In ^ we compare our results with 
one-boson-exchange potential evaluated using the nucleon-nucleon-meson coupling obtained 
in our previous paper .'3' Section [H] is for a brief summary. In the appendices, we review the 
ADHM construction of two instantons and summarize the necessary formulas used in this 
paper. 

§2. Nuclear force in holographic QCD 

In this section, we briefly summarize the treatment of the single baryon in the D4-D8 
modeP''^ of the holographic QCD following Ref. flSjl and describe our strategy for obtaining 
the nuclear force. Our first goal is to obtain a quantum mechanics Hamiltonian for a two- 
nucleon system. The total Hamiltonian consists of one-body canonical kinetic terms for 
each nucleon, potential terms for each nucleon, plus interactions. One generically has an 
interaction potential as well as a correction to the kinetic term. Then, secondly, we evaluate 
the inter-baryon energy using the Hamiltonian. This provides an explicit nuclear force that 
is dependent on nucleon states labeled by spin and isospin. 

The concrete calculations of the Hamiltonian will be given in ^ and its evaluation with 
explicit nucleon states will be presented in detail in ^ 

2.1. Review: single baryon in the model 

First, we review briefly the single baryon case^*' in the holographic QCD proposed in 
Refs. ED, [in])- The notation of our paper follows that of Ref. [T5]) . 

Our starting point is the meson effective action derived in Refs. [9|l. [T0|) . which is given by 
the following five-dimensional U{Nf) Yang-Mills-Chern-Simons (YMCS) theory in a curved 
background: 

S = Sym + Sqs } 

Sym = —k, I d!^xdz tr 



Scs = ^, I u:M) ■ (2-1) 



247r2 

Here /x, z/ = 0, 1, 2, 3 are four-dimensional Lorentz indices, and z is the coordinate of the fifth 
dimension. The field strength is defined as J-" = | J-'a/jfix" A dx'^ = dA + iA A A with the 
U{Nf) gauge field A = Aadx"" = A^dx'^ + Azdz (a = 0, 1, 2, 3, z), and the front factor k is 
related to the 't Hooft coupling A and the number of colors as 
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The action (12- ip is written in the unit Mkk = 1, where Mkk is the only dimensionful 
parameter in the modelO The functions h(z) and k{z) appearing as the "metric" in the 
action (12- II) are given hj h(z) = {1 + z'^)~^^^ and k{z) = 1 + z"^, while, in the second term, 
u^i^A) is the CS 5-form (here, we omit the A product, e.g., AJ-''^ = A A A J-") 

u,{A) = ti(^A^'-'-A'J'-^A'^ . (2-3) 

In the two-flavor case {Nj = 2) that we focus on in this paper, the U{2) gauge fields A are 
decomposed as 

A = A + A— = A''— + A— = y" A^— , (24) 
2 2 2 ^2 ^ ^ 

c=o 

where t"- [a = 1, 2, 3) are Pauli matrices and r° = I2 is a unit matrix of size 2. 

This action is obtained from the low-energy effective action on Nf D8-branes in the 
curved ten-dimensional geometry corresponding to Nc D4-branes wrapped on a circle with 
an antiperiodic boundary condition for fermions. At low energy, this D-brane configuration 
provides U{Nc) QCD with Nf massless quarks and the action (12 -ip describes the dynamics 
of mesons and baryons. The action (12- ip is written in (1+4) dimensions, and the space 
along the extra dimension x'^{= z) is curved. Once the gauge fields are decomposed into 
their Kaluza-Klein states concerning the z direction, each mass eigenstate corresponds to 
a meson, and the action (12- ip describes the whole spectra/interactions of the mesons. By 
contrast, baryons are solitons with nonzero instanton number in the four- dimensional space 
parameterized by = [x, z) (M = 1, 2, 3, z). As they are localized in the four- dimensional 
space in the five-dimensional space-time, they behave as pointlike particles. The instanton 
number is identified with the baryon number® and these particles are interpreted as baryons. 
We will see more of the details below. 

The single-baryon solution was found to have the size of order X-^/'^]iSinB jg helpful 
to rescale the coordinates as^^' 

~M _ \l/2 M ~0 _ 

Aoit, x) = Ao{t, x) , Auit, x) = X-^^^Auit, x) , (2-5) 

to see the consistent 1/ A expansion of the equations of motion and the total energy of the 
single baryon. Hereafter, we omit the tilde for simplicity. Then, for large A, the YM part of 



*' In Refs. iSj) andHni), these two parameters are chosen as Mkk = 949 MeV, k = 0.00745 to fit the 
experimental values of the p meson mass rup ~ 776 MeV and the pion decay constant /tt — 92.4 MeV. 



6 



the action is 



'S'ym = — clNc I d'^xdz tr 
^ ' d'^xdz 



A ^2 
2 ^MN 



— 
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while the total equations of motion are 



(2-6) 



^mnpqFmnFpq + 0{\ ) — 



D^Fmtv + 0{\-^) = . 

dnFoM + 77;; — ^^MNPQ S i^iFMNFpq) + —FmnFpq 

dr,FMN + OiX-') = . 



(2-7) 
(2-8) 
(2-9) 

(2-10) 



Therefore, at the leading order, the warp factors h{z) and k{z) are approximated by 1, so 
the SU{2) part of the equations is nothing but the standard YM equation in flat space. It 
is solved by a BPST instanton located around z ~ 0. The electric U{1) part is sourced by 
the instanton density, as seen in fl2-9p . The explicit solution is^^-^ 



ACI 



^0 



with the BPST instanton profile 



1 - 



2^2 



Aa = A 



M 







(2-11) 



e 



Z)+i{x-X) ■ T 



z-zy 



If -XI 



(2-12) 



p is the size of the instanton, while = (X^,X^,X^, Z) = (X, Z) is the position of the 
soliton in the four- dimensional space. 

Quantization of this soliton has been carried out in Ref. [TSll . It is basically the same as 
the quantization of a YM instanton in the moduli space approximation,'^''^ except for the 
additional potential in the moduli space induced by the presence of the subleading terms 



in the action fl2-6p . The moduli space for a single YM instanton is Ai^ 



X M7Z2 



parameterized by {X,Z) and (/ = 1,2,3,4) with the Z2 action — )■ —y^. The radial 
component p = a/ {y^Y of y^ gives the instanton size and the angular components = y^ / p 
parameterize the SU{2) orientation of the instanton. The quantization of the soliton is 
described by quantum mechanics on this moduli space, with the Lagrangian 



'fx^ + ^Z^ + ^iyr-Uip,Z) 



(2-13) 
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where the "mass" for each moduh is given by 



mx = mz = ^ = STi^aN, , (2-14) 



and the potential 



Uip. Z) = M„ + 8 AiV. + ^1 + |!) (2.15) 

is obtained by substituting the solution (12-111) to the action (12 -Op . Mq = Stt^k is the classical 
mass at the leading order in the 1/ A expansion. The potential is classically minimized at 

which shows that in fact the soliton has the size of order A~^/^ when it is rescaled back to 
the original coordinates by (12- 5p . The Hamiltonian is given by 

-1 / d V -I f d V -I f d ^"^ 



This system has an S'0(4) ^ {SU{2)j x SU(2)A/'^2 rotational symmetry acting on . Here 
SU{2)j and SU{2)j are interpreted as the isospin and spin rotations, respectively, and they 
act on y = + iy'^T^ as 

y^giygj (2-i8) 

with {gi,gA ^ SU{2)j x SU{2)j. The isospin and spin operators are given by 

2 dy- y dy^ """'^ dy-) ' 

respectively. From this, we have P' = and, hence, only baryons with I = J appear in this 
approach. 

Quantum states of the baryon can be labeled using quantum numbers of isospin/spin 
I = J = 1/2 , (/ = 1, 3, 5, ■ ■ ■ ), the eigenvalues of the third components of isospin and spin 
operators and J^, and the quantum numbers Up = 0, 1, 2, ■ ■ ■ and = 0, 1, 2, ■ • ■ , which 
label the excitation numbers of (almost) harmonic oscillators in p and Z, respectively. For 
example, the proton and neutron have quantum numbers {I, l3,np,nz) = (1,1/2,0,0) and 
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(/, Ja, np, Uz) = (1,-1/2, 0, 0), respectively. The corresponding wavefunctions are normalized 
spin/isospin stated^ 

1 i 

\p t) = -iy' + w^)/p , b i) = — (y^ - w^)/p , 

TT TT 

\n t) = -{y' + W')/P , \ni) = --{y' - W')/P , (2-20) 

TT TT 

multiplied by the following p and Z wavefunctions, 

R{p)=p'e-^^\ MZ) = e-^'\ (2-21) 

with [ = -1 + 2^1 + Up = 1/VQ, and uz = ^J^Jl. The functions R{p) and ^z[Z) 

should be multiplied by normalization factors. 

2.2. Our strategy 

Our strategy for the calculation of the nuclear force consists of three steps: 

1) Construction of generic two-baryon solution of the YMCS theory ( 12-ip . 

2) Computation of the quantum-mechanical Hamiltonian for the moduli parameters, and 

3) Evaluation of the Hamiltonian with specified nucleon states. 

This is a direct generalization of the single-baryon case to the two-baryon case. In the 
following, we describe each step in more detail. 

2.2.1. Construction of two-baryon solution 

The case of two baryons, which is our concern, can be considered by solving the equations 
of motion of the original action (12- ip with the constraint that the instanton number is 2. As 
we have seen, the rescaled variables are useful for seeing the properties of the single baryon. 
There, one can start with a BPST instanton solution in flat space, since the size of the 
instanton is smaller than the scale of the curved background geometry. When we have two 
baryons, the situation is different. If the distance between the two is larger than O(1/Mkk)0 
(or 0{\/X/Mxk) in the rescaled coordinate), the effect of the curved space-time comes into 
play, thus a similar analysis cannot be performed. In this paper, we concentrate on the case 
where the two baryons are close to each other, i.e., the distance is smaller than (9(1/Mkk)- 

It is well-known that one can explicitly construct generic two-instanton solutions of Eu- 
clidean four-dimensional YM theory in flat space. We use ADHM construction of the in- 
stantons for our purpose. The construction is reviewed in Appendix [Bl 

The two-instanton moduli space is parameterized by four quaternionic parameters {Xi, X2, 
y\i 2/2)- We summarize our notation for the quaternion in Appendix |X1 The quaternion has 

As a reference, if we use the value of A/kk that is fixed by the rho meson mass, we have 1/Mkk — 
0.208 fm. 
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a representation by 2 x 2 complex matrices as in (lA-31) . In this notation, these moduh 
parameters can be written as 

Xi = Zi + tXi-T , y^ = yf + iyi-T, (2 = 1,2) (2-22) 

with Xi = {Xl,Xf,Xf) and -jji = {yl , yf , yf) . When the separation between the two instan- 
tons is large, the two-instanton solution can be approximated with a superposition of two one- 
instanton configurations with moduli parameters {Xi,yi) {i = 1,2). Here, X^'^ = {Xi,Zi) 
corresponds to the position of the instanton in the four- dimensional space, pi = a/ y( y- is 
the size, and = fji/pi is the SU{2) orientation of the instanton. 

Defining r*^ = Xf ^— X^^^ (M = 1, 2, 3, z) as the relative position of the two instantons and 
|r| = V r^r^' as the distance between them in the four-dimensional space, the requirement 
for the fiat space approximation to be valid amounts to 

|r| < C (^v^/Mkk) . (2-23) 

Note that this is written in the rescaled coordinates fl2-5p . 

As seen from the structure of the equations of motion in the 1/A expansion, the only 
nonzero quantities at leading order are Am and Aq, as in the case of the single baryon. The 
equation of motion fl2-9p shows that the U{1) part of the gauge field is again sourced by 
the instanton density, now with two maxima at the location of the separated baryons. The 
explicit solution of the SU{2) two-instanton solution and the U{1) part will be presented in 
§31 with the help of the ADHM construction reviewed in Appendix |Bl 

2.2.2. Calculation of quantum-mechanical Hamiltonian for two-baryon moduli 

The next task is to obtain the classical potential U{y{, y2, X1 — X2, Zi, Z2). We substitute 
the two-instanton configuration into the action f l2-6p . As we mentioned, the nonzero fields 
at the leading order are only Am{x) (the spatial components of the SU{2)) and Aq{x) 
(the temporal component of the U{1)). Therefore, in the rescaled action ( l2-6p . nonzero 
contributions are 



U = 2Mo + Hltt'^'^^ + hIT'^ + O (A-i) , (2-24) 

J d^xdz tr [z^FIjj^] , (2-25) 



i/^ff ^^^^ = aNr I d^xdz tr 



f2 



(2-26) 



The first term of the total potential f l2-24p is the energy contribution from the first term in 
the action fl2-6p . which gives the leading order term in the 1/A expansion. The important 
part is the subleading order potential, H^^l^^'^^^ and H^l;^^\ which, in the case of the single 
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baryon, was used to determine the size of the instanton classically and was responsible for the 
quantum effect. We compute these for the cases with two baryons, with explicit dependence 
on the moduli parameters. In the last equation of (12-251) . we have used the self-dual equation 
for the instanton. 

The computation of Hp^^^"^^^ is straightforward although it is technically involved, which 
will be presented in §3.11 On the other hand, it turns out that H^j^^^^ is not easy to compute. 
One can evaluate it numerically, but numerical results are not useful for our purpose, as we 
will later need to make a moduli integration of it with the baryon wavefunction. Thus, we 
need the explicit analytic expression for the moduli dependence of the potential. For this 
purpose, we concentrate on the case with a large inter-baryon distance, 

0(1/Mkk) < |r| , (2-27) 

where the left-hand side is the size of the single instanton, fl2-16p . in the rescaled coordinate. 
In this region, since there is only a slight overlap of the instantons, one can obtain an analytic 
expression for H^l;^^\ The evaluation will be presented in §3.21 

Together with the constraint fl2-23p . in this paper, we consider the separation of the 
baryon satisfying 

O (1/Mkk) < |r| < 0{VX/Mkk) (2-28) 

in the rescaled coordinates fl2-5l) . 

The quantum mechanics of the moduli parameters consists of the potential term U and 
the kinetic term. The kinetic term of the quantum mechanics of the moduli is given by the 
moduli space metric. As opposed to the single-instanton case, the moduli space metric of 
the two-instanton configuration is complicated. It is found that the asymptotic form in the 
case of a large separation takes the form 

ds'^ = dsl + dsl + 0{\r\-^) , (2-29) 

where dsj = 2{dy{Y + (dXfY + 2{dylY + {dX^^)"^ is just two copies of the metric for the 
single instanton, and dsl is (9(|r|^^), which is our concern. This contributes to the quantum 
mechanics as a (9(|r |^^) correction to the kinetic term of the moduli dynamics. We explicitly 
compute this correction in §3.31 

All together, the analytic expressions of the (9(|r|^^) terms of the quantum mechanics are 
computed, and this is the Hamiltonian for the two-baryon interaction. The total expression 
is summarized in §3.4[ 
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2.2.3. Evaluation of nucleon-nucleon potential 

The final step is to evaluate the energy with the Hamiltonian of the quantum mechanics. 
The interaction Hamiltonian is given in 1/1^1 expansion, and we treat this as a perturbation. 
The state of our baryons is specified at infinite separation, and we use simply the tensor 
product of two copies of a single-baryon wavefunction. This can be justified at leading order 
in the perturbation of quantum mechanics, *3 

Since we have an analytic expression for the Hamiltonian, the integration of the moduli 
parameter with the given wavefunctions is straightforward. The result is the nucleon-nucleon 
potential, in particular if we choose the baryon wavefunctions to be that of a nucleon. This in- 
tegration will be presented in §11 The distance between the baryons in the three-dimensional 
space, |r| =a/ (XJ— X2)^ + (Xf — + — Xf)^, is related to the four-dimensional distance 
\r\ as \r\ = ^/\r\'^ + {Zi — Z^)^. Note that we fix one of the moduli \r\ and perform the inte- 
gration of the other moduli, Zi, Z2, ?/(, y\. This is because we are interested in the potential 
in the scattering problem, rather than the computation of the bound state energy. 

Our final result for the nucleon-nucleon potential is given in fl4-45p and (14 -46 p . The 
central force f l4-45p shows that the nucleons have a repulsive core. We also obtain the tensor 
force f l4-46p . All the potentials have the form which is peculiar to four- dimensional 

space, as described in the introduction. 

To illustrate the properties of our nuclear force fl4-45p and f l4-46p . we next compute 
the one-boson-exchange potential among nucleons. In our previous paper,'^'' we derived the 
nucleon-nucleon-meson coupling in the D4-D8 model of the holographic QCD. By using this 
coupling, the summing up of all types of mesons propagating among the nucleons should 
provide a certain aspect of the nuclear force. This computation can be carried out for 
arbitrary distances between the nucleons, as long as the nucleon radii do not overlap each 
other. The resultant nuclear force, at larger distances, exhibits the standard properties of 
the nuclear force, such as scalar/tensor forces due to pion/p- meson/other-meson exchanges. 
The computation will be presented in ^ 

We will find there that, in the region (12 -28 p . this one-boson-exchange potential does 
not coincide with our nuclear force (I4-45P and (I4-46P derived using the ADHM construction 
of two instantons. The reason is that when nucleons are close to each other the nucleon 
itself is deformed by the effect of the other nucleon. In deriving the one-meson-exchange 
potential, this effect cannot be taken into account. Thus, naive computation based on the 
one-boson exchange is not sufficient to capture the complete picture of the nuclear force at 

*' For describing the deuteron system, this perturbation is not the way to proceed. One needs a minimum 
of the whole potential of moduli including r, to obtain quantized energy of a bound state of two baryons. 
Our interest in this paper is the potential force appearing in the scattering process of the two baryons. 
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short distances. 



§3. Effective Hamiltonian for two baryons 

In this section, we calculate the effective Hamiltonian for the quantum mechanics of the 
two-baryon state following the strategy described in the previous section. The system is 
described as a quantum mechanics of a particle living in the two-instanton moduli space. 

3.1. Potential from SU{2) part 

Let us first evaluate the contribution of the SU (2) part of the gauge field to the potential. 
As explained in ^ the leading term in the 1/A expansion can be obtained by substituting 
the two-instanton solution in the fiat space-time into the action fl2-6p . The two-instanton 
solution can be obtained using the ADHM construction. See Appendices |A] and |B] for our 
notation and a brief review of the ADHM construction. 

The leading term in the SU (2) part is obtained by evaluating (12 -25 p . This integral can be 
calculated using the formula flC-7p obtained in Appendix O Substituting flB-14p into flC-7p . 
we obtain 

j d^xdz tr Fli^ = Svr' {pj + pl + 2{Zf + Z|) + 4:{w^Y) , (3-1) 
where w"^ is the real part of (IB-17p . which is given by 



w 



' - ' {yivt - mvt + my<yi) = ^^tr {tr'^a.'a,) . (3-2) 



|^|2 ^^-^^^ ■ »^ aij 2 |^|2 

Therefore the potential from the SU{2) part is given by 
KT''^ = ^ {pI + pI + 2{Zl + ZD + Ain^^f) 



pi + pi + 2{Z't + Zl) + plpf-^ ii{iT^a^^ai) ii{iT'a^''a^) ) . (3-3) 



r 



The leading r-independent terms reproduce the contribution of SU{2) part in the one- 
instanton potential fl2-15p for the two instantons. The next-to-leading term gives the in- 
teraction between the two baryons. It is of order 1/|tP as expected in a five- dimensional 
gauge theory. 

3.2. Potential from U{1) part 

The field strength of the SU (2) part of the gauge field satisfies (IC-3P and then it is easy 
to see that 

^o = 7^niogdetL (3-4) 
327r^a 
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is the regular solution of the equation of motion (12-91) . which vanishes at infinity. Here, L is 
given by (]B-2p and (1B-I4p . By using the constraint (1B-I6p . it can be written as 



L{x) 



fi{x) e{x) 
e(a;) /2(x) 



where 



fi{x) = p1 + \x-Xi 



\w\ 



1,2) 



e(x) = {yi ■ y2) + {w ■ (Xi + X2 - 2x)) . 



(3-5) 



(3-6) 
(3-7) 



We evaluate the energy contribution of this U{1) gauge field when the separation between 
the two instantons is large. Then, it can be confirmed that the energy density is mainly 
concentrated around a; ~ (i = 1, 2), where the two instantons are located. When x is 
close to Xi, we can make the expansion 



X ~ Xi, \x — Xi \ <^ \x — X2I . 



(3-8) 



Without losing generality, we can choose Xi = and X2 to be very far away from the origin, 
and so we choose x around the origin. The order estimate is 



1X0 



I X — X2 1 ~ I r I = I X2 — Xi I ^ I — Xi I ~ I I 



(3-9) 



We can evaluate the matrix L in this expansion, and obtain the following expression: 

4 



□ log det L 



\x\ 



1 



Pi 



+ 



1 + 



{\x\^ + piy 

2{yi ■ y2?pl 



{\x\^ + plf 



\x? + plf 



\w\^ + 0{\r\ 



In this expression, note that w = 0{\r\ ^). More explicitly, from (]B-17p . we have 



T^m X yi\ 



(3-10) 



(3-ii; 



Then, the gauge field around the position of one of the two instantons x ~ Xi is expanded 
as 



A 







\x\ 



1 



p\ 



{W + p\? 



1 



1 



2Ypl 



{W + pIY 



+ 0{\r\ 



(3-12) 



where 



Y = {yi- y2f - \y2 x = 2(2/1 • - \yi?\y2\^ = pIpI (2(ai ■ "2)' - l) • 



(3-13) 
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The leading term in the '\-/\r \ expansion reproduces the result (12- lip for the single instanton. 
We are interested in the next-to- leading term of order \r\~'^. Note that the first term in the 
next-to-leading terms, 

1 



1 



(3-14) 



is precisely the leading contribution of the five-dimensional Coulomb interaction with the 
second instanton located at a; ~ X2. 

Let us compute the potential energy, with this expression for the gauge field. The U{1) 
part of the energy (12-261) can be written as 



pot 



(fxdzAQUAQ . 



(3-15) 



Note that to obtain this expression we performed a partial integration, it is harmless at 
this stage because the gauge field Aq decays fast asymptotically. Then, we substitute the 
expanded expression f l3-12p to this energy formula. This procedure is slightly ambiguous, 
since the integration does not commute with the l/\r\ expansion in (13-120 . We give a more 
systematic way of evaluating the integral in Appendix [El Here, we present an easy way to 
obtain the correct answer. Substituting the expansion (I3-I2p into (I3-I5p . we obtain 



407r2a 



Pi 



Pi 



(3-16) 



as the leading term. This is just the sum of the energy contribution of the single instanton 
in (12-151) . The next-to- leading order that we are interested in is the cross terms, 

2Ypl 



2 

aNr, 



d xdz- 



d xdz 



1 

(87r2a)2 



\x\ 



Pi 



1 + 



\xr + Pi 
2Ypl 



2\2 



{\x\^ + plf 



□ 



□ 



1 + 



(1=^1 



\x\ 



p\ 



\x\ 



(3-17) 



Here, (7/1 O 7/2) denotes the contribution from the integration around x ~ X2, which is 
obtained by exchanging y\ and y2 in the first and second terms of (I3-I7p . Note that the 
first term in (I3-17P is different from the second term. This is because, at this stage, the 
partial integration suffers from a surface term, owing to the constant l/|Tp. Anyway, we 
can perform the integration analytically, and the result is 



pot 



A^. 



407r2a 



1 



1 



2 2 

p{ pi 



+ 



1 2{a^-a2? 

2 5 



1 



2 2 

PI pi 



+ 0(|r 



-3\ 



(3-18) 
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3.3. Kinetic term 

The kinetic term of the quantum mechanics of the two-baryon states is given by 



mx 



(3-19) 



where X° = (Xf ,?//) 

space that are promoted to time- dependent variables and X'^ 
the two-instanton moduh space. 

The hne element of the two-instanton moduli space for large separation is obtained in 
Ref . BSD as 



[a 



1,2,..., 16) are the coordinates of the two-instanton moduli 

4X". is the metric of 



ds^ = dsl + dsl + 0{\r\ 



(3-20) 



where 



dsl 



dsl 



{dXi ■ dXi) + {dX2 ■ dX2) + 2{dyi ■ dyi) + 2[dy2 ■ dyi) 
2 r 

plidyi ■ dyi) + pl{dy2 ■ dy2) + 2{yi ■ dyi){y2 ■ dy2] 



- (?/2 ■ dyif - {yi ■ dy2f - 2{yi ■ y2){dyi ■ dy2) 
+ 2eijKLy{yidy^dy2 - ((2/2 ■ dyi) - {yi ■ rft/s))' 



(3-21) 



The leading terms ds^ in (13-211) is just a sum of the metric for each instanton, which gives 
the canonical kinetic term f l2-13p . The next-to-leading terms dsf contribute to the potential 
of order l/|rp. 

The kinetic term of the Hamiltonian is given by the Laplacian of the two-instanton 
moduli space as 

1 _o 



kin 



2m 



X 



The outline of the calculation is summarized in Appendix The result is 

= V^ + V? + 0(|r|-=^) , 



(3-22) 



(3-23) 



where 



d 



dXf^ 
1 



+ 



d 



9Xf 



+ 



Pi 



d 

dyi 



2 yOyi 



d 

dyi 
2 



+ 



y2 



d 

dyi 
2 



dyi 



yi 
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(3-24) 



d 



d d 



+ eijKL yivig^g^ + {vivi + y^i - ivi ■ y2) 5'') 



dyi 

d d 
dyi dyi 



(3-25) 
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Using the spin/isospin operators fl2-19p . it can also be written as 



d 



M 



d 



d 



2 \pldpi y^dp 



pI 



i=l,2 
4 



1 d 



PI dpi V ' dpi 



d 



Pi 



P: 



ja ja 

2 i i 



(3-26) 



p\ 



A Tarn d d d d f J d 

dpi dp2 dpi dp2 V 9y{ 



Vi 



(3-27) 



Here, a,b,c = 1, 2, 3 are the indices for the SU{2) adjoint representation and the subscripts 
i = 1,2 label the two instantons. 

3.4. Summary 

By collecting ([S3]), f E^ . and flH^ with flH^ or (jMID, the total Hamiltonian is 
obtained as 



H = Ho + Hi. 



(3-28) 



The leading order Hamiltonian Hq is just two copies of the Hamiltonian for one baryon fl2-17p 
obtained in Ref. USD 



i=l,2 



-1 



d 



+ 



1 / d 



-1 



d 



U{pi,Zi 



2mx \dXiJ ^ruy \dyjj 2mz \dZ, 

where U{pi, Zi) is the potential given in fl2-15p . 

The final term Hi gives the (9(|r|^^) interaction between the two baryons 

1 



(3-29) 



Hi = iff + H[ 



r{SU{2)) 



2m 



X 



where mx = Sn'^aNc and 



Nc 1 



1 2{ai . - 1 (pI p\ 

P\ pI 



(51/(2)) 



and V? is defined in dMZD- 



^— ^PiP2-r4j tr(ir"a2 ^Oi) tr(ir^a2 ^Oi^ 



(3-30) 

(3-31) 
(3-32) 



§4. Nucleon-nucleon interaction 



We are ready for the evaluation of the nucleon-nucleon interaction potential at short dis- 
tances, using the quantum mechanical Hamiltonian (I3-28I) obtained in the previous section. 
This section is devoted to demonstrating these manipulations in detail. 
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We evaluate the interaction Hamiltonian Hi with the nucleon wave function obtained in 
Ref. [T5]l . The wavefunctions for consistent quantum states of two asymptotic nucleons are 
given in §4.H and with the wavefunctions, the computation of the expectation value of Hi 
follows in §4.21 Finally, the decomposition of (Hi) into a central force and a tensor force is 
given in §4.31 Our final result for the nucleon-nucleon potential at short distances is fl4-45p 
and fl¥iB]) . 

4.1. Wavefunctions of two-nucleon states 

We have two baryons, and the wavefunction for them is given by a tensor product of the 
wavefunctions fl2-20p and fl2-2ip . We can arrange f l2-20p in a simple form, 

1. 2 ^ / bt) bi) . .... 



IJ 



where is the Pauli matrix with a = 2. Here, we specify the matrix element using 
the indices I, J that take values in {±1/2}. From (14- ip . we see that {I, J) component 
of the matrix r^a is directly related to the wavefunction of the spin/isospin state with 
{P, J^) = (J, J). For the two nucleons specified by {If, Jf, /|, J|), our wavefunction is 



\iT^ai)i,j,{T''a2)i,j, , (4-2) 



where we have omitted the upper index 3 in If and Jf for simplicity. 

The wavefunction for the instanton size variable p is given by a tensor product of (12-211) . 
which is R{pi)R{p2). To reduce the computational effort for the evaluation of the p integral, 
we use the following simplified wavefunction instead of (12-211) . 

R{pi)R{p2) , where R{p) = p' exp [- p^] . (4-3) 

Hereafter, we mean R{p) using this expression. This is obtained by a rescaling {l/2)myUjpp'^ — )■ 
p^. Note that, among the terms in Hi (13-301) . ijj^''^^'' and are invariant under the rescal- 
ing, so we can just replace the wavefunction by this simplified one. However, ff^^'^^'^^^ has 
a scaling dimension that needs to be taken into account. The overall normalization of (14- 3p 
will be taken care of later. 

The wavefunction for the Z modulus, for the nucleon states of = 0, is again a tensor 
product of (l2-2ip and given by 

^lJziZl)MZ2) - (4-4) 
This is, again, up to a normalization constant. 
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4.2. Spin/isospin matrix elements of interaction Hamiltonian 

Let us evaluate the expectation value of the interaction Hamiltonian Hi, for given wave- 
functions for bra and ket, (14- 2p . (14- 3p . and fl4-4p . Since the moduli Zi and Z2 are decoupled 
from the rest, we treat them later independently. Here, first, we integrate the moduli y{ and 
?/2- The integration measure is d}yid'^y2 = plpldpidp2df]idf]2, where dili is the integration 
over the angular coordinates af in SU{2) ~ S^. The normalization is given by J dfli = 2-7?^, 
which is the volume unit S*^. For the integral df2i, the following formulas for the integration 
over g G SU(2) are useful: 

f -1 -1 

J dn g.,g,,'gmn9,,' = -[2 {5u5^,5,,5^p + 5.,5^,5M - 5a5^,5,,5^u - 5,u5np5.,5^i] , 
j dQ gijg^l = n'^SuSjk ■ (4-5) 

Evaluation of h[^^^^'^ 

Let us evaluate the expectation value of h['^^^^^ f l3-3ip . In h[^^^^\ the spatial coordinates 
X^^ do not couple to y-, so it provides only a central force, and does not yield a tensor force. 
As for the integration over the p variables, only the following three integrals appear: 

/•OO /"OO 

dpp'Rip)\ / c^pp^i^(p)^ / dppR{pf. (4-6) 
Jo Jo 

On the other hand, to compute the expectation values, we need to include the normal- 
ization of the wavef unctions. Considering both, we only need the following ratios for our 
computation: 

J^^dp p'Rip)'^ I [J^dp p'Ripy^ = 1 + ^ ' 

(^J^^dp pR{pr^ I [J^dp p'Ripf^ = . (4-7) 

These can be derived by partial integrations. Using these, we obtain 



/ dpidp2 p\pI f 4 + 41 R{pifR{p2) 
J \Pi P2/ 



p!^p!/ jdpidp2 plplR{pifR{p2Y "^FTT' ^^'^^ 

2 2 

Note that we could use f l4-3p because ^ + ^ is scale-invariant. 

Next, let us perform the a integration. In Hf^^^^\ the nontrivial term is only (ai ■ 02)^, 
so we compute the matrix element of this. We first note 

(ai • a2f = ^ tr a\a2 tr aia\ = ^{a\) k 1(0.1) mn (0.2) LK{al)NM ■ (4-9) 
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Thus, with the wavefunction (14-21) . the matrix element is 
((ai • a2f)^j,^j^j^^j^)^^^j^j^^j^) 



Here, we have used ((r^a)// 



I2J2 



(4-10) 



a'T For the integral f l4-10l) . we can use the formula 



fl4-5l) to obtain 

(4-11) 



To obtain the last expression of fl4-lip . we have used the identity 

Combining fl4-lip with fl4-8p . we obtain the matrix element for H^^'^^'^^ as 



(4-12) 



HI 



{u(i)) 



Nc 1 
Bvr^a Irl 



(/(,j;,/^,4),(/i,Ji,/2,J2) 
1 3/ + 1 



10 [+1 



11 + 2 



(4-13) 



The Zi, Z2 dependence in |rp will be integrated later, together with the other terms in Hi. 
Evaluation of ^^^^^ 

Next, we consider fl3-32p . First we make the integration over pi and p2- In the present 
case, the integral is not invariant under the scaling of p as opposed to the previous examples. 
This time, after the scaling {l/2)myUJpp'^ — > p^, we obtain the multiplicative new factor 
A/ {niyUpY . Therefore, using (14- 7p . we obtain 



{pIpD 



i /rfpi p\R{pif jdp2 plR{p2) 

J dpi p\R{piY jdp2 plR{p2] 



Next, we consider the integration over ai and 02 • It is obvious that it proceeds in the 
same manner as use the formulas f l4-5p . we bring the relevant part of uf^^"^^^ to 

the form 

tr(ir"a2 ^ai) tr(ir^a2 ^"1) = tr(r"a2 ^oi) tr(r^a];^a2) 

= {ai)ij{a:^^)pq{a2)qr{a2^)ki{T°')jk{T^)rp . (4-15) 



1 + 2 



(4-14) 
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Then, using (14-51) . we obtain 



= g [26ab{4:6i'^i^6i;^i.,6j'^j,6j;^j2+6i'^i.,6i^i^6j'^j^6j^j2) 
To simplify this expression, we use (14 •121) and the Fierz transformation 



(4-16) 



(4-17) 



Then, fl4-16l) becomes 



(tr(ir''a2 ^ai) tr(ir^a2 



Jl,h,J2) 



J'J2 



+ g(^/(/i ■ ^/^/2) ((^'')4J2(^^)j(Jl + (^'')j^J2(^")j(Ji-^ab(rj^J2 ■ Tj(Ji)) 

Combining (14-141) and (14-181) . we arrive at 

jj(SU{2))\ 

^ I (/(,J(,/^,4),(/l,Jl,/2,J2) 

(^ + 2)^ |rT r, , , c 



(4-18) 



(4-19) 



where r = rj\r\. 
Evaluation of 

Finally, we evaluate in Hx- In (13-271) . only the last two terms {v\'§^ + (^29^) 
have the angular dependence. Thus, as for the terms other than these, we only have to 
perform the p integral. For example, using 



1 A (p^^R^p) 

dp \ dp 



r. 1 



4p^ + (-8-4/) + (/^ + 2/)-)i?(p) 



(4-20) 
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the ratio formulas (14- 7p can be applied, and this results in 



,Ji,h,J2) 



(i+2){!+5) 



l + l 



.1 d 



(4-21) 



(I[,J[,I^,J!,),{Il,Jl,l2,J2) 



We have used 4/"/f = f/j/^ • t/^/j- 

To perform the integral in the last term in ( l4-2ip . we first perform a partial integration 



.1 d 



V'(/l,Jl,/2,J2) 



Using the wavefunctions (14-21) and (14- 3p . we find 



dyi 



(4-22) 



I d I I 



(1-^02)/! Ji (r^a2)/2 J2 —R{pi)R{p2) 
Pi 

dpi V Pi 



^(P2 



(4-23) 



Thus, in fl4-22p . the following integrals for p appear: 

jpidp^jpidpARipirRip^r 1+2 



I Pfdpi J pldp2 R{pifR{p2) 
hldpij Pldp2 plR{p2?R{pi)^^ 



l + l 

d_ ( Rjpi) 
pi 



jpldp,jpidp2 R{piyR{p2Y 

I Pldpij pldp2 plplR{p2f 



a 


{"¥)] 


dpi 





j p\dpij pldp2 R{piyR{p2f 



(-2)'^, 

l+l 

il + 2){!+5) 
l + l 



(4-24) 
(4-25) 
(4-26) 



Here, again, we have used (14 -Tp . Then, straightforward calculations with the formulas (14-5 
and fl¥TTi) show that 



.1 9 



(/ + 2)(/ + 5) 



dyi 



(I[,J[,I^,4),{h,Jl,l2,J2) 
2^I[h^I^l2^JiJl^4J2 + Y^^^^'l^l ' ^^2^2)(^J(Jl ■ ^Jp2' 



(4-27) 
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As a result, the matrix elements of the third term in fl3-30p are obtained from fl4-2ip and 
fl¥?n) as 



2m 



X 



yi)ii[,j[,i^,4Uh,Ji, 



,h,J2) 



327r2aA/'J rP 



(1 + 2) (1 + 5) 



1 + 



+ 5) f 1 \ 



(4-28) 



Since / ~ 0{Nc), this term is 0{1) in the 1/Nc expansion, while the other terms f l4-13p and 
(14- 191) are 0{Nc). Therefore, this term is subleading, compared with the contributions from 
iff and //p^(^)\ 

Summary of the result 



Summing up all the results fl4-13p . fl4-19p . and (14 -28 p . we have the spin/isospin matrix 
elements of the interaction Hamiltonian 



1 \r\'^ 

(^i)(/(,j(,7^,j^),(/i,ji,72,j2) = ^1 + ^2 



(4-29) 



where the coefficients are given by 



Cl 



87r2a 



+ 



10 / + 1 



45 / + 1 



327r^aN, 

il + 2){l + 5) 



l + l 



^i[h^i'^i2^j[.hh!,j2 - a^^i[h ■ ^/^/2)(^j;ji ■ ^Jp2) 



C2 



([+2)2 
327r2aA^^ 



5i[hh'^i2h'^j^h^J2 + g(T7(/i • rr^h) (2(r ■ T)j^jAr- f)jjj, - fj^j^ ■ tj'^j^ 



(4-30) 



(4-31) 



4.3. Final result for nuclear force at short distance 

We finally evaluate this Hamiltonian fl4-29l) with the wavefunction for Zi, (I4-4I) . and take 
the leading term in the large Nc expansion. 

In our result (l4-29p . the Z dependence is included in the four- dimensional distance, 
|r| = a/I?^^ + {Zi — ^2)2. We fix the baryon distance |r| in the real space, and perform the 
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integration over Z\ and Z2. To perform the integration, we rewrite the wavefunction f l4-4p 
as 



exp 



(4-32) 



Since ( I4-29P consists of two terms, l/|Tp and l/l?"!^, we use the following formula for inte 
gration, 

7re"'^'(l -Err[a/3]) 



00 ^-fp-x"^ 



dx 
dx 



e-- - 20Fa/3 + 7r(l - 2a'^(3^)e'^''^\l - Err[a/3]) 



Here Err is the error function. 



Err \x\ 



2a3 



e"*' dt . 



Then we obtain 



v/?/3e^'l^1'/^(l -Err[/3|r1/2]) 
2^^ ' 
/3|r1 + v/^(l - /3^|rlV2)e^'l"1'/^(l - Err[/3|r1/2]) 

4|^3 



/3 



(4-33) 
(4-34) 

(4-35) 

(4-36) 
(4-37) 



where /3 = ^2mzujz = (16(2/3)i/2^2^Ar,)i/2. 

For large Nc, the argument of the error function, /3|r| is very large for nonzero |r|. We 
can use the following asymptotic formula for the error function. 



1 _ 2 
1 - Err(x) = —=e 



X Ix-^ 



-4\ 



By using this expression, the Z-integral results are markedly simplified, 



1 



(4-38) 



(4-39) 



z V? \klVz 1^' 

This expression can be easily guessed since this is just a substitution of Z\ = Z2 = 0. The 
value Zi = Z2 = is the classical value of the location of the instanton in the z space, and 
the large limit should reproduce the classical result. 

On the other hand, the coefficients ci and C2 in fl4-29p are evaluated in the large Nc 
expansion as 



(81 3 \ 

— '5/(/i^/^/2^J(Ji<^4J2 + -(T/^/i ■ Ti^I2){tj[J, ■ Tj^J^)] + 0(1) , 



(4-40) 
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+2(77,,, ■ f,.,J(f ■ r)j(j,(f ■ t)j>j,) + 0{l) . (4-41) 

Here, for deriving these, we used = {A/5)N^ + 0{1). It is interesting that in fact all the 
contributions from Vf drop off, since they are subleading compared with the terms from 
H^^^^^'^ and nf^^'^^'^ in this large N^. expansion. 

Therefore, the leading term of the interaction Hamiltonian is 



-l^5i[h5i'^i2^j[.h5j'^J2 + ^(^/(/i • T'i'2i2W- r)jf^j^{f- t)j>j^ j . (4-42) 



(-^l)(/(,J(,/^,J^),(/l,Jl,/2,J2) 

_ vriVc /27. . . . 6 
~ W 

We rescaled r back to the original coordinate (remember the rescaling f l2-5p ) . This is the 
nucleon-nucleon potential at a short distance in the large limit. 

Let us decompose this force into a central force Vc(|r|) and a tensor force Vrd^*!), 

V = Vcm + SuVrm . (4-43) 

Here, the tensor operator 5*12 is defined by 

Su = 3(ai • f )(a2 ■ f ) - • (72 = 12(Ji ■ f )(J2 ■ h - 4/1 ■ h , (4-44) 

where Oi = {a}, af, af) = 2 Jj are the Pauli-spin operators (it is just twice the spin operator). 
Applying the decomposition (14-431) to our result (14-421) . we obtain 

^c(|r1) = vr f ^ + ?(™(«)1 , (4-45) 



^T(|r1) = y/r/2y ^ - (4-46) 

We have derived the nuclear force at a short distance. As already mentioned, the ex- 
pression is valid in the region (12-281) . We have found that there is a strong repulsive core in 
the central force (14-451) . Our finding is quite important, as an analytic computation of the 
repulsive core of the nuclear force, based on strongly coupled QCD. 

The tensor force (I4-46P is found to be negative for / = (since /"/2 = —3/4), which is 
consistent with the experimental observation and also with lattice QCD calculations.'^ It is 
intriguing that the nuclear force has |r|~^ dependence. This \r\~^ is peculiar to the physics 
with one extra spatial dimension, and thus it is a direct consequence of the holographic 
approach to QCD. It would be interesting to fit the lattice and experimental data of the 
nuclear force with our result. 
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§5. Comparison with one-boson-exchange potential 



In this section, for a comparison, we compute the nucleon-nucleon potential using the 
one-boson-exchange picture. This can be done by integrating the meson propagator with 
the nucleon-nucleon-meson coupling obtained in our previous papei® in the D4-D8 model in 
the holographic QCDO We find some discrepancy between the two pictures. The reason for 
this discrepancy is basically the fact that in the one-boson-exchange picture the deformation 
of the nucleon by the other nucleon has not been taken into account, while our potential 
obtained in ^ includes this effect via the ADHM construction in the previous section. In 
sum, we find that the one-boson-exchange potential captures merely a part of the nucleon- 
nucleon potential. 

5.1. Interaction potential 

In Ref.EI), the nucleon-nucleon-meson couplings were computed, by extracting the asymp- 



totic behavior of the one-baryon solutionol The one-boson-exchange potential is obtained 
by just evaluating the energy of a superposition of two asymptotic baryons, which are re- 
garded as point particles sourcing the meson fields propagating between them. The distance 
r should be larger than the instanton size p, of course, because we use the asymptotic form 
of the solution (which was given in §2.3 in Ref. |S|) as a baryon solution. 

Before getting into the details of the evaluation of the potential energy of our system, it is 
instructive to consider a simple example of a two-electron system in classical electrodynamics 
with the action 

S = J d'x(^-\F'^^-fA,^ , (5-1) 

where is a current of an external source. Here, we work in the Lorenz gauge dfj,A^ = 
and consider a static configuration. Then, the equation of motion AA^ = is solved by 

A^{x) = {A-^fm ^ [ d'yA-\x,y)riy) , (5-2) 



where A ^{x,y) = -^j^zrg^ is the Green's function for the Laplacian A. The on-shell action 
is evaluated as 

S = -\l d'xA^AA, = -^l d'x^A, = -^l d'xrA-'j, . (5-3) 

In ^only, we use the upper (or lower) index "(1)" and "(2)" to label the gauge fields and the currents 
for the two instantons, and we do not use the rescaled coordinates (|2-5p . This is for making use of the 
notation of Ref. ^ . 

**^ In Refs. the couplings were computed by a different approach, using baryon spinor fields in 

the bulk curved space-time. 
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Let us consider the current associated with two electrons placed at Xj {i = 1,2): 



j«(x)=e5=^(f-X0, jf =jf =jf = 0. (. = 1,2) 
Then, the solution of the equation of motion is given by 



(54) 



l(2) 



A 



(5-5) 



Substituting this in the on-shell action fl5-3p and picking up the cross term, we obtain the 
potential due to the interaction of the two sources as 



(5-6) 



Let us follow this line of argument for our action f l2-ip in the curved space-time. It was 
shown in Ref. [8) that nonlinear terms in the equations of motion can be neglected in the 
asymptotic region. Therefore, we are allowed to consider the linearized equations of motion 
in the curved space-time: 



hiz)dlAi + d,ikiz)d,Ai) = , 9jA + d,ih{z)~'d,ik{z)A,)) = 
with the gauge condition 

h{z)d^A^ + d,{k{z)A,) = . 
The potential energy analogous to (15 -Op is 



V = 2k (fxdz tr 



(2) 

where the "current" ja is defined as 



(5-7) 



(5- 



(5-9) 



j^;^^ = {h{z)dA + dMz)d.)A^o^ 



jf = {h{z)dA + dMzmA) 



(2) 



(5-10) 



Here, Aa^ {i = 1, 2) are the asymptotic solutions for a single baryon located at x'^^ ~ 
that satisfy the linearized equations of motion (15- 7p and the gauge condition (l5-8p . They are 
exphcitly obtained in Ref. E]). The "current" (15-lOp behaves as the pointlike source placed 
at Xg^, which corresponds to the delta function source jj^'^ in (15 4 p in the previous example. 
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Of course, the asymptotic solutions can only be trusted in the asymptotic region, since the 
nonlinear terms in the equations of motion will become important near the position of the 
instanton X^^ . However, here, we simply assume that the baryons can be treated as point 
particles to obtain the one-boson-exchange potential. 

For simplicity, we consider the static configuration, *3 By noting that our focus is on the 
leading order behavior of the large A and large N^. limit, the leading contribution to the 
currents turns out to come only from 

lo,^,^, (5-11) 

which are of order 0{\^^) while the next-to-leading order is by the component Aq, which 
is of order 0{\~^N^'^). All the other components are of order 0{\~'^N~^). Thus, let us 
consider only the leading order components. 
The solutions and the "current" are 

= ^^'i^ - - Z,) , (5-12) 

-2Ti\p,)hT {T'a,r\a,)-') [f""'^ - 6'^^^ 5\x - X,)5{z - Z,) , (5-13) 



,■(2)6 

Ji 



= -2n\p,)\i (rVr"(ai)-^) ^Hix, z; X,, Z,) , 

= -2n\p,)hT {r'a.T'^ia.r') ^S'{x - X,)6{z - Z,) . (544) 

Note that as in Ref . , these expressions are written in terms of the original variables without 
the rescaling f l2-5p . keeping them in the order of 0{1). Here, G and H are the Green's 
functions associated with the linearized equations of motion ( \5-7\\ . They are obtained in 
Ref. ED as 

oo oo 

G = Kj2MZ2)MZi)Ynm) , H = f,Y,MZ2)MZi)Ynm , (5-15) 

n=l n=0 

where {i'n}n=i,2,- is a complete set of the eigenfunctions satisfying the eigenequation 

- h{zy^d^{k{z)d^i)n) = Ktpn , (5-16) 



This amounts to throwing away momentum-dependent nuclear force such as L ■ S force. 
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with eigenvalue A„ and the normahzation condition 

K j dzh{z)i)ni)m = 5nm , (S'lT) 

and {4>n}n=o,i,- is defined as 

(Pn{z) = -^d,iJn{z) , 0o(^) = u . ■ (5-18) 

The eigenfunction ipn{z) is an even (odd) function for odd (even) n. When the five- 
dimensional gauge field is expanded by {ipniz)}, the coefficient fields are interpreted as 
the vector (for odd n) and axial-vector (for even n) meson fields, and the eigenvalues A„ are 
proportional to the mass squared of the corresponding mesons. The function y„(|r|) is the 
Yukawa potential associated with the vector/axial- vector mesons of mass \/X^ : 



1 p-V>^\r\ 

Ynm) = -7 ^ • (5-19) 

Att \r\ 

We substitute these expressions to the interaction potential (15- 91) to obtain 
1 



V =K 



-G{X2, Z2] Xi, Zi 



4a2A2 

+ An\p^f{p2)hI (rVr"(ai)-^) tr (r'a^r^a^)-') 

^ (^^"^'^ - '^w^ i^'^'^i - '^w) 

+ 47r^(pi)2(p2)'tr [r'a,r\a,)-') tr [r' a^r^a^Y') x ^^//(X^, Z^; Xi, Z^) 

(5-20) 

This is the inter-baryon potential energy that we want to evaluate as the nuclear force, in 
the one-boson-exchange approximation. 

5.2. Short distance behavior 

When the distance between the solitons is smaller than I/Mkk, the Green's functions 
G and H can be approximated by their fiat-space analogue as explained in our previous 
paper,'^ 

G = H=^^ ^ i . (5-21) 

47r2|Xi-X2P + (Zi-Z2)2 

Substituting this expression to the inter- instanton potential energy fl5-20p . we can easily see 
that only the first term in (l5-2Up remains, while the other terms cancel each other and vanish, 

V = . (5-22) 

IGttVA^ |Xi - X2I2 + (Zi - ^2)2 
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This is a harmonic potential in four-dimensional space. Once the 2;-part of the wavefunction 
for the nucleon, fl4-4p . is taken into account, the expectation value is given just by substituting 
the classical value Zi = Z2 = to the above potential f l5-22p . 

^ = — f^- (8-23) 

Here, |r| is the distance between the nucleons. We find that only the central force appears. 

This expression (15-231) is apparently different from our previous results (14 -4511 and (14-461) . 
In fact, in (l5-23p . there is no isospin dependence, and no tensor force. The reason is that the 
one-boson-exchange description is not sufficient to capture the entire nuclear force. In fact, 
the two-instanton solution in ADHM construction is not just a superposition of two BPST 
instantons. In particular, there is a deformation of the instanton that causes additional 
contribution. The one-boson-exchange potential (I5-23P only captures a part of the complete 
results (I4-45P and (I4-46P obtained in the previous subsection. 

In Appendix [Fl we try to evaluate the potential height for nucleons overlapping each 
other in real space, in the one-boson-exchange picture. 

5.3. Large distance behavior 

Let us look at the large distance behavior of the nucleon-nucleon potential obtained from 
(l5-20p . The large distance means \Xi — X2I ^ 1 in the unit Mkk = 1- In this limit, 
essentially only the pion dominates, since only the pion is the zero mode while others have 
Yukawa potential that decays exponentially fast. 

The contribution of the pion comes from n = component of the Green's function H 
defined in (I5-I5p . Therefore, in the potential (l5-20p . we are interested in the last term, and 
the function H is now approximated by a massless Green's function in three dimensions: 

V~ K- 47r'{pi)'{p,)hi (T»a,T"(oi)-') tr (tVt'Io,)-') 

For spin 1/2 baryons, the trace part can be easily evaluated as 

64 

tr (r^air"(ai)-i) tr [r' a2T%a2)-') = -{lll^Jtr^ , (5-25) 

y 

where Jj and /j are spin and isospin operators, respectively. The potential energy is then 
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where Su is defined in f l4-44p . Here, we have used the relation 



d d 



3 

1^ 



(5-27) 



with f = Xi — X2 andr = f/\r\. 

This expression can be compared with the well-known one-pion-exchange potential 



1 / fl-A 



m 



S^2 ( m. 



4 



(5-2^ 



In the chiral limit m^r — )■ 0, only the tensor force remains, and it agrees with fl5-26p when 

(5-29) 



\k{Z)/ 



U 3 

which is exactly the relation found in Ref. 18]). 

If we use the classical values for the above expectation values 



/ (pi? \ I ip2) 



\A;(Zi)/(i) \k{Z2 
then the potential (15-261) becomes 



— Pel 



(2) 



16 



Svr^aA V 5 ' 



(5-30) 



V 



2Nr 



d d 



9Xf 9X||Xi-X2 



■(A%')^i2^. (5-31) 



307r2aA 



In a quantum evaluation for the expectation values, we will have roughly x (1.05)^ times the 
classical value above, after substituting the numerical valuesP 

5.4. Intermediate distance behavior 

As the baryons approach each other from asymptotics, there appear effects of the massive 
meson exchange. At this intermediate distance, the potential (15-201) becomes 



4a2A 



^ 00 



n=l,odd 



25fi7r'* 

(J f) (j (j 



+ 



2567r^ 



dXi dXl dZ^ dZ2 



n=l 

DO 

J2 MZ2)MZl)Yr,{\f\) 



71=0, even 



(5-32) 
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where we have used fl5-15p and also dropped the terms including {d/dX){d/dZ), tp2k{Z), 
and (j)2k-i{Z) with k = 1, 2, ■ ■ ■ , because they vanish when they are evaluated with baryon 
wavefunction, owing to the parity property of the ipn{Z) functions. The potential can be 
summarized to the following form: 



V ^ K 



+ 



4a2A2 
2567r4 



9 



J2 MZ2)MZl)Yn{\f\) 



71=1, odd 



X 



9 

5 d 

dXfdXi 



J2 MZ2)MZi)Y„i\r]) + S'"'J2 KMZ2)MZi)Ynm 



(5-33) 

n=0,even n=2,even 

The first term in the first line, l/{4a'^X'^), corresponds to the contribution of the vector mesons 
that appear from the U{1) part of the gauge field. Among them, the u meson exchange is at 
the lowest order {n = 1). The second term in the first line is the contribution of the vector 
meson in the SU{2) part of the gauge field, whose lowest order {n = 1) corresponds to the 
p meson. The third line gives the contribution of the pion {n = 0) and axial-vector mesons 
[n > 1), whose lowest order is the ai meson, in the SU{2) part of the gauge field. Note 
that the contributions from the U{1) part of the axial-vector and pseudo-scalar mesons are 
sub leading in the 1/Nc expansion. 

To divide this expression into the central force and tensor force, we use the following 
formula for the Yukawa potential y„(|f^). 



_d d_ 

dX? dXi 



■Yum =[—- ) ( ^ + ^ + An ) K(|rl) - — A„r„(|r1) . (5-34) 



3 



3 



1, 



Then, we obtain the potential energy due to the central force and tensor force, 

V=Vc + S^2Vt , 
1 



(5-35) 



Vc=k' 



+ 



oo 



4a2A2 
2567r^ 



^„(^2)v^n(^i)v;(|ri; 
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n=l,odd 

{p,nP2)\iiii){JiJ2) 



2 

X - I - 
3 



A„V^„(Z2)^„(Zi)y„(|r-1) + Y A„0„(Z2)0„(Zi)y„(|rX 



n=l,odd 



n=2,even 



(5-36) 



,647r' 



iPi)\p2f{llll) ( E ^n{Z,)MZi) + ^ + A„) Y^{\r\) 

\n=l,odd ^' ' ' ' ^ 



n=0,even 



3v^A^ 
|r| 



+ \n]Y^{\r\] 



(5-37) 
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We note that these expressions can be reproduced from the leading order of the one-boson- 
exchange potential computed from tree-level Feynman diagram with the nucleon-nucleon- 
meson couplings obtained in Ref. |8]). See Appendix O for detail. 

Let us look at the contribution from light mesons. The u meson gives only the central 
force, the first line in Vc with n = 1. It is 

Vt'^ = -«:^^(^i(^i))(i)(^i(^2))(2)ri(|r1) . (5-38) 

For the potential between the same types of baryons, this expression of course reduces to 

V!,^^ = -^'^,{UZ)?Y,{\r\) . (5-39) 

Since the Yukawa potential y„ defined in fl5-19p is negative, the u meson exchange gives a 
strong repulsion force in the central force. The u meson is the lightest vector meson that 
comes from the U{1) part of the five-dimensional gauge field. The instanton is electrically 
charged under this U{1), so the baryons should have this universal repulsive force. The u 
meson exchange manifests its lowest term in the KK decomposition of the higher- dimensional 
"electric" repulsion. 

The p meson exchange can be seen in SU{2) components of n = 1. The central force is 
^^,) = «:^^(pi)^(p2)^(«)(JiV2") MZ2)MZi){-Xi)Y^m , (540) 

while the tensor force is 

Kp) _ .26471^ _ f 3 , 3v^ 



When two baryons are of the same type, these reduce to 



= ^'^{pr{nn) {uz)y ^^+^+^1)^1(1^1) • (5-41) 

The strength of this tensor force can be compared with the strength of the pion tensor 
force (15-261) . The ratio of the front coefficients is given by —{ipi^Z))"^ / (0o(^))^- The classical 
evaluation of this gives 

'u'i^y^I = -^^(MZ))' ^ -7r(0.597)2 ~ -1 . (5-42) 

Let us see how this p meson exchange may change the result of the vr exchange. At the length 
scale r ~ 1 fm ~ 200 MeV"\ the pion behaves as a massless particle while the p meson is 
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massive. The above ratio denotes that the coefficient is of the same order, while the sign is 
opposite to each other. The normahzation of the Yukawa potential gives the following rough 
ratio 

f,-mp\r\ 

- ^(0-1) • (5-43) 

This means that the p meson exchange does not contribute much to the tensor force at the 
distance scale ~ 1 fm. 

§6. Summary 

In this paper, we have deduced the nuclear force at short distance in large A^^^ strongly 
coupled QCD, by applying gauge/string duality. 

In the D4-D8 modeP''^ of holographic QCD, baryons are instantons in (l+4)-dimensional 
YMCS theory. We have explicitly constructed a two-instanton solution in the theory by em- 
ploying ADHM construction of instantons. The analytic expression for the potential energy 
plus kinetic terms of the baryon, i.e., the effective Hamiltonian of quantum mechanics for 
two-baryon system, has been derived, for the distance C(1/(V^Mkk)) <r < C(1/Mkk)- 

The evaluation of this interaction Hamiltonian for specific two-nucleon states, labeled by 
spin Jf and isospin If with i = 1,2, provides the nuclear force at the short distance scale. 
We have obtained a central force f l4-45p as well as a tensor force f l4-46p . The central force 
exhibits a strong repulsive core of a nucleon. As the repulsive core has been mysterious from 
the viewpoint of strongly coupled QCD, our result is of importance as a derivation of the 
repulsive core from the "ffist principle" of QCD, in the large Nc expansion and also with the 
gauge/string duality. 

The obtained nucleon-nucleon potential at short distances has behavior. Technically 
speaking, this comes from the harmonic potential in four spatial dimensions including the 
holographic extra dimension. It would be interesting to fit this peculiar behavior with the 
experimental observation or the recent lattice result .'^ 

As our result is for the short distance, it is important to generalize our analysis to a larger 
distance scale. For r > (9(1/Mkk), the effect of the curvature along z becomes indispensable, 
thus, one needs to construct a two-instanton solution in curved space. This may lead to an 
analysis of a deuteron system in holographic QCD. For this, the inclusion of quark mass to 
the modeP^ that should make the pion massive may be important. On the other hand, the 
height of the nucleon-nucleon potential at r = is of interest, but its derivation has turned 
out to be difficult, as we explained in §2.21 Further effort along these directions may reveal 
some more interesting physics in QCD, via the holography. 
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Appendix A 

Notation for quaternion 

A quaternion q G EI is given as a linear combination of the form 

q = q'^-q^I -q^J -q^K , (g" G M , m = l,2,3,4) (A-1) 
where J, J, and K satisfy 

= = = -1 ^ IJ= -JI = K , JK= -KJ = / , KI = -IK = J . 

(A2) 

Since (— ir^, — ir^, — zr^) satisfy the same algebra as (/, J,K), a. quaternion can be represented 
as a 2 X 2 complex matrix of the form 

q = q'^ + iq^T" . (A-3) 

The conjugate and norm of a quaternion q are defined as = -|- q^I + q^J + q^K and 
\q\ = a/ q'^q = \fq^ = y/q'^q"^ , respectively. 

The product of two quaternions q,w E M follows from the relation (1A-2P and it can be 
written in terms of the 2x2 complex matrix representation flA-31) 

qw = q^w^ — q- w + i{q^w + w'^q — qxw)-T, (A-4) 

where q = (g^, g^, g^), etc. We also use the following notation: 

{q-w) = ^{q^w + w^q) = q"'w"' , (A-5) 
{q X w) = -{q'^w — w'^q) = i{q^w — w'^q + qx w) ■ f . (A-6) 
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An element of Sp{n) is defined as a n x n quaternionic matrix Q satisfying Q'^Q = ln *J 
In particular, an element of Sp{l) is a quaternion satisfying 

q^q = q^q^ = 1 , (A-T) 

which is equivalent to the condition for an element of SU{2) in the 2x2 complex matrix 
representation (lA-3p . 

Appendix B 

ADHM construction 

B.l. ADHM construction for Sp{n) instantons 

Here, we briefly review the ADHM construction."^ (See for example Ref. [32]) for a 
review.) Since SU{2) = Sp{l), the ADHM construction for the Sp{n) instantons is useful 
for our purpose. 

We define an [n + k) x k quaternionic matrix of the form 

A{x) = a + bx, (B-1) 

where a and b are {n + k) x k quaternionic matrices and ^ J — x'^K G H is a 

quaternion composed of the coordinate of the four-dimensional space (x^, x^, x^, x^) = {x, z). 
The matrices a and h are required to satisfy that a^a and 6^6 are kxk real symmetric matrices, 
and a^h is a k x k symmetric quaternionic matrix. These conditions for the matrices a and 
b are equivalent to the constraint that the matrix A satisfies 

A^A = L{x) (B-2) 

with akxk real symmetric matrix L{x). We also implicitly assume that the matrices a and 
h are generic and they are matrices of rank k. 

The ADHM gauge field for the /c-instanton configuration is given by 

Am{x) = -iU{x)^dmU{x) , (B-3) 

where U{x) is {n + k) x n quaternionic matrix satisfying 

Z\tf/ = , f/tf/ = 1„ . (B-4) 

Note that the matrix U (x) is defined up to a transformation 

U{x) -> U{x)g{x) , {g{x) G Sp{n)) , (B-5) 
*' Sp{n) in this paper is the unitary symplectic group, which is also written as USp{2n). 
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which acts as a gauge transformation for the gauge field (lB-3p . 

The gauge field (IB-Sp as well as the constraint (]B-2|) is invariant under 



A{x) QA{x)R , L{x) R' L{x)R , (B-6) 

where Q G Sp{n + k) and R G GL{k, R). By using this transformation, the matrix b can be 
fixed as 

b=\ ° 1 , (B-7) 



-1 



k 



and then A is of the canonical form 



k 

where X is a k x k symmetric quaternionic matrix and Y is an n x k quaternionic matrix 
such that Y'^Y + X'^X is a k x k symmetric real matrix. Note that we have not completely 
used the transformation (IB-6p . In fact, a transformation (IB-6P with R G 0{k) and 

Q = I ^ I , e Spin)) (B-9) 



leaves flB-7p invariant. 



B.2. = SU(2) one-instanton 

As an exercise, let us consider the n = k = 1 case. Using the canonical form (IB-Sp . we 
have 

(B.10) 

with y,x, X E H. In this case, (IB-2P is satisfied without imposing further constraints. 
The condition (IB-4P is solved by 



where ^ = y^\x — and p = 



Then the ADHM gauge field flB^ is 
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where 



x-X 



9 



(B-13) 



and a = y/p is an element of Spil) = SU{2). 

This is the BPST instanton solution. The moduli parameters X, p, and a correspond to 
the position, size, and gauge orientation of the instanton. The a-dependence of the gauge 
field can be eliminated by a global gauge transformation. However, it is known that this 
degree of freedom is also physically relevant when we quantize the system via moduli space 
approximation method. 

B.3. Sp{l) = SU{2) two-instanton 

For n = 1 and k = 2, the ansatz flB-8|) can be written as 



Y 



[yi,y2) 



X 





w \ 


, A{x) = 




\ w 


X2 J 





yi 

Xi — X 
w 



y2 \ 

w 

X2-X J 



The constraint (1B-2|) requires 

Y^Y + X^X = 
to be a real symmetric matrix and hence 



\w\ 



\yi\' + \Xi\- 

ylvi + Xlw + w^Xi 




y\y2 - ylyi + r^w - w^r 







where we have defined r = Xi — X2. This equation is solved when 



w 



[y2 xyi) + ar 



(B-14) 



(B-15) 



(B-16) 



(B-17) 



with a G M. This parameter a can be eliminated by the residual 0(2) symmetry in ( \B-9\i 
with q = 1 and R G 0(2)^^ We will set a = in this paper. 

After all, we have 4 quaternionic parameters yi, y2, Xi, and X2 to parameterize the two 
instanton moduli space. It can be shown that when the separation of the two instantons is 
sufficiently large, the two-instanton configuration can be approximated by the superposition 
of two 1-instanton configurations. Then, Xi {i = 1,2) represents the positions of the two 
instantons, pi = a/Ij/iI^ corresponds to their size, and = yi/ Pi is their SU{2) orientation. 
This fact can be explicitly seen in the effective Hamiltonian (13 -28 p . 
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Appendix C 

ASL 
'pot 



Evaluation of H^^^^"^^^ 



In this section, we compute 

jd'xzHiFl^, (C-1) 
where z = x^. We will follow the strategy of Ref. in which 

j d'x\x\HiFl^, (C-2) 

is calculated. 

To evaluate trF^^, we use the useful formulaP'' 



tr Fl^ = ^e™ tr FmnF„ = -□□ log det L , (C-3) 



where □ = dmdm and L{x) is defined in (lB-2p . 
In general, L{x) can be written as 

L{x) = A\x\^ - 2^mx'^ + A , (C-4) 

where A = b'^b and A = a^a are positive definite k x k real symmetric matrices, and 7^ 
(m = 1,2,3,4) are k x k real symmetric matrices. For the canonical form (IB-Sp . they are 
given by 

yl = l, 7m = ^m, A = Y^Y + X^X, (C-5) 

where is the k x k real symmetric matrix satisfying X = X4 — Xil — X2J — X^K. 
The result we are going to prove is 



d'xzHiFl^ = SttHt [{i,A-'y - (71/1-^)' - {i2A-'y - (73^-^)^ + AA^'] . (C-6) 
For the canonical form with f lC-5p . we obtain 

j d^xzHiFl^ = S'R^ii{2{X^f + Y^Y) . (C-7) 
Note that this formula implies 

j d^x\x\^tiFl^ = lQ'K^ii{2Y^Y + X^X) , (C-8) 
which reproduces the result in Ref. |3i 
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To show (]C-6|) . we can set yl = 1 without loss of generahty. The A dependence can easily 
be recovered by the transformation (lB-6p . Integrating by parts, we obtain 



d X z ti 



— lim 



d-'QR'x 



2 



z'^dm.O log det L — 2z5m4: □ log det L + 2dm log det L 



(C-9) 



where = and d^f2 is the volume element of unit S^. Here, we have used Gauss's law 



d'^x dm W„ 



lim 

-R->-oo 



(C-IO) 



S3 



where Um = x"^/R is a unit normal vector of the 5''^. 
We can easily check the following formulas 

9m log det L = tr{Vm) , 
□ logdetL = 8tr(L-^) -tr(V^) , 
a^DlogdetL = -12tT{L-^Vm) + 2tT{VmK 



n I ' 



(C-11) 
(C-12) 
(C-13) 



where Kn = ^dmL. Using these, flC-9p becomes 



d X z ti 



— lim 



d^QR^ tr 
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- Uz^L'^x^^Vm + 2z^ x'^VmV^ 



IQz'^L-^ + 2z^V^ + 2x"'Vrr, 



(C-14) 



Inserting the relation 

d^xzHiF^^ 

d^QR^ tr 



L-^dmL = 2L 



7m), we obtain 



lim 



53 



-1 



(C-15) 



+ (L + 2{R^ - x™7m)) {R^L~^ + L-^lmL-^lm - 2L--^x^^m) 

We are only interested in the 0{R^) terms in the integrand of the right-hand side of this 
equation. Then, recalling L = 0{R^), we obtain 



d^xz^ tiFL 



lim 

R-^00 



d^Q tr 



S3 



Pi + P2 + Pz + Pa + P5 



(0-16) 



where 



Pi=R^ ■ AL-\R^ - Az^) , P2 = ■ 4x"7„ L~\-5R^ + W 



P3 = R 



16(x"7n)'(2/?' - 3z^) 



R- 



8A{R^ - Qz^ 



P^ = R-^ ■ 2A-iy 



(C-17) 
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Useful formulas to evaluate the integral are 



53 



TC 



S3 



S3 12 



(no sum for m, n) 



(C-18) 



To evaluate the integral of Pi in (lC-17p . we expand L 



-3 



(C-19) 



and 



lim / d-'ntiPi 
^53 



- hm / £n tr 

J53 

- lim / d^Q tr 

J53 



+ 4(i?^-4z^)(-AA+^(x-7„.r 



(C-20) 



Using the formulas (IC-lSp . we see that the first term vanishes. The second term also vanishes 
since it is odd under x — )■ —x. From the third term, we obtain 



lim / d^f2 tr Pi = -16n^ tT{-fl + -f^ + 7I - 87, 



(C-21) 



The integrals for the other P2, ■ ■ ■ , P5 are calculated in a similar manner The results are 



- lim / d^n tr P2 = Svr^ tr (4(72 + 7^ + ^l) - 372) , 

- hm f ci3^trP3 = -47r2tr(3(7? + 72' + 7|)+74) , 

-lim / d^f2 tTP4 = 7r^tT{A) , 
J53 

- lim / d^f] trPs = -127r2tr(7^) . 

J53 

Summing up all these, we finally obtain 

I d'x tr Fl^ = Svr^ tr (-(7? + 72' + l!) + ll + A) . 



(C-22) 

(C-23) 
(C-24) 
(C-25) 



(C-26) 
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Appendix D 

Laplacian of the two-instanton moduli space 



Here, we outline the derivation of the expression (13-251) and fl3-27p . 
Consider a metric given as 



(D-l) 



where g^i^ is a constant metric and h^p ^ 1 is a small perturbation. Then, omitting the 
0{h'^) terms, we have 



(D-2) 



1 



(D-3) 



where (^"^) is the inverse matrix of g^/^, Vg = 'g'^^dadp, g = det(5'Q,^), g = det(^^^). Here, 
raising and lowering the indices is done with ^"'^ and We apply these formulas to the 
metric f l3-2ip . in which ds^ and dsi correspond to ^^.^ and hap, respectively. 
A little algebra shows 

2 



{d,h\)d^ 



\r\ 
2 



2(l?/il + \y2 

d 

yi 

Vi 



dyi 
d 

dyi 



+ I 2/2 

2/2 



d 

dy2 

d 

dy2 



Oi\rn 
0(\r\~^) 



(D-4) 
(D.5) 

(D-6) 



d d 



pI 

2 \dyi dyi 
dyi 



+ 



d d 



- 2/2 



pI 



2 \ dy2 dy2 

2 



2/1 ■ 



_d_ 
dy2 



- (2/1 • 2/2) 



2/1 



d 



dyi 
d d 



2/2 



_d_ 
dy2 



dyi dy2 



From these, we can easily obtain fl3-25l) . 

The following formulas are useful to obtain the expression in (13-271) : 



2/i ■ 



_d_ 
dyi 



_d_ 

'dpi 



(D-7) 



(D.8) 
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d_ d_ _ 

dyi dyi) dpj ' pidp. 



•9' 3 9 4 
+ —- -If If 



•9' 3 9 4 

+ -T^ ^JfJf 



dpi Pi dpi P: 



2 « J ' 
i 



ja ja 



ja ja 



{yi ■ y2) 
ivi ■ y2) 



d d 



dyi dy2 
d d 
dyi dy2 



.i„.j 



d d 



d d 



(D-g) 

(D-IO) 
(D-ll) 



Appendix E 

EvaluaUon of HlT^ 



Here, we rederive the result fl3-18p in a more systematic way. From the expression in 
fl3-5p . we obtain 



where 



log det L = log /i + log /2 + log /s 



Mx) = pi + \x-X,\''+\w\\ {t = l,2) 



f3{x) 



h{x)f2{x) 



with 



e(x) = (yi ■ ya) + (w ■ (Xi + X2 - 2x)) . 
Substituting (13 ■41) into fl2-26p . we obtain 



mi)) 

pot 



2 (327r2a)2 
aNr 1 



d^xdz {dMD{\ogh + log/2 + log/3))' 



d xdz 



J2 (dMa log + {dMD log 



i=l,2 



2 (327r2a)2 

+ 2(9Mn log /i) (^Af □ log /2) + 2 ^ (Sa/D log /,) (Sa/D log /a) 

1=1,2 

The following formulas are useful for evaluating this integral: 
□log/, 
□□log/, 
□□□log/. 



A{\x-X,\ 


' + 2(p? + 1' 


u;|2)) 


(|a;-X,|2 + p2 + 
96(p2 + |iy|2) 


w 

2 


2^ 


{\x-X, 


? + p1 + 


w 


2-)4 ' 



1536(p2 + 


w 


2)2(3 


X- Xi 


2 




m + 


w 


'))) 


{\x 




' + p1 + 


w 


2))6 



(E-l) 

(E-2) 
(E-3) 

(E-4) 



(E-5) 



(E-6) 
(E-7) 
(E-S) 
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for i = 1,2. 

Then, the first term in flE-Sp is evaluated as 



J d^xdz {dMU\oghf = - j d^xdz {U\ogh){UU\ogh 



1 r r 2 ,4(^2 + 2) 96 



pI + \w\^ Jo (m2 + 1)2 (m2 + 1)4 

1 2567r2 



pI + \w\'^ 5 
2567r2 / 1 \w\^ 



Pi Pi 



0{\rn , (E-9) 



where we have used 



u^^==l=, (E-10) 



and M = lii-l. A similar formula for /2 is obtained by replacing pi with p2 in (IE-9p . 
The third term in (]E-5|) is evaluated as 



rf^xdz(9Mniog/i)(9Mniog/2) = -y rf^xd;2(niog/2)(nniog/i) 

1 



4(|\/iw + f|2 + 2lf ) 


96 


{\Viu + + v^y 


+ 1)4 



where r = Xi — X2 and we have defined 



Vp1 + 




\r\ 





f=^, ■ (^ = l,2) (E-12) 

To evaluate the leading term in the l/\r\ expansion, we consider the limit — )■ 0. Although 
the integrand of flE-lip is divergent a.t u = —r/Vi when V2 = 0, the integral around u = 
—r/Vi is convergent. Besides, there is a suppression factor 1/(m^ + 1)^ that makes the 
contribution around u = —r/Vi in the integral vanish in the Vi — )■ limit. Therefore, we 
can safely take the Vi limit and using 



/ 



4 ■ 96 

d\j^^^^=64.\ (E.13) 



we obtain 



/ 647r2 
d'xdz{dMDlogh){dMa\ogh) = + 0{\r\-') . (E-14) 
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The last term in (lE-Sp is given by 

d'xdz (9Mniog/i)(9Mniog/3) 

d\dz {UUU\ogh)\ogh 



PI + IIOP J + 1) 

with 



1536(3^^ - 2) , 
^ " — / 9 , log /s , (E-15) 



((2/1 • 2/2) - 2Ki 


r (lu 


u)f 


r 




f Viu 


2 + 



/3 = 1 - ^...r ^.2^ • (E-16) 



Note that we have used the relation {w ■ r) =0, which follows from the definition (1B-17P 



with a = 0. Again, to obtain the leading order terms in the 0{\r\ ^) expansion, it is allowed 
to pick up the leading term in the integrand as 



Using the formula 



,4 1536(3^2- 2) 1 1287r2 ^ 



(m2 + 1)6 ^2 + I 



we obtain 



/ d'xdz{dMn\ogh){dMn\ogh)^^^j^^^{a^-a,y + 0{\r\-') , (E-19) 
J 5 \r\^ PI 

and similarly 

/ rf=^a:rfz(aMniog/2)(aMniog/3)^^^^4(«i-"2)' + 0(|rr=^) . (E-20) 
J ^ |r| P2 



As one can see in flE-17p . log/3 is (9(|r| and hence the second term in flE-51) does not 



contribute to the leading 0{\r\ terms in the potential. Collecting ( ]E-9I) . (]E-14p . (]E-19p . 
and (lE-20p . we reproduce the potential fl3-18l) . 

Appendix F 

Height of one-boson-exchange potential 



In this appendix, we try to evaluate the height of the nucleon-nucleon potential, in the 
one-boson-exchange approximation. Note that as shown in §5.31 the one-boson-exchange 
model does not describe correctly the short distance behavior. Thus, this appendix is only 
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for an illustration of what will happen in general when two baryons are on top of each other 
in real space. 

If the instantons are located within the distance of (9(1/Mkk), the one-boson-exchange 
potential is f l5-22p . When the instantons are located at Zi ^ 0, there is an additional classical 
potential coming from the self-energy part fl2-15p . so in total, the inter-instanton potential 
energy is 



V = , — ^ + STT^aXNc 



(F-l) 



This classical potential exhibits an interesting structure. Let us find a minimum of this 
potential for fixed inter-baryon distance in real space, |Xi — X2I = |r|. We employ a classical 
approximation for Zj, by taking the large Nc limit, for simplicity. Owing to the exchange 
symmetry Z\ ^ Z2, the potential is minimized at Zi = —Z2 = ^4/2. Thus, the minimization 
problem is for the potential 

V = : M^. , .2 + —^^l ■ (F-2) 



IGvr^aA |?^^ -|- r| 



The minimization condition is 



dV r4 Sn^aXN, 



dr4 Svr^aA (|r|^ + r]^ 

This is solved with 



This forms a three-dimensional sphere around the origin in the four-dimensional space. For 
a fixed |r|, we obtain nonzero r4 to minimize the classical potential. The instantons go away 
from Z = axis, to minimize the potential energy. This can be understood as follows. The 
instantons have the overall U{1) electric charge, so they try to be away from each other. 
But at the same time there is an effect of the curved space-time, which tries to bring the 
instanton toward the Z = axis. The balance of these two effects results in the minimization 
at f lF-4p . The minimum energy for a fixed inter-baryon distance is 

^3 3 ' ' ^ ^ 

Thus, the inter-baryon potential height is maximized at \ f\ =0, with the height value 

Knax = ^ . (F-6) 
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The sphere (1F-4|) does not reach the region |r|^ > In fact, in this region, the 



potential energy is minimized by r4 = in flF-2p . which results in the previous result fl5-23p . 
It is smoothly connected with fIF-Sp at |r|^ = g^^ . 

For the sphere (lF-4p to make sense, its radius should be larger than the classical radius 
of the instanton, (12- 161) . Unfortunately, both are of the same order, ~ 0{1/ vStt^vA), so one 
cannot trust this sphere radius. However, we find an interesting picture of the instantons, 
where the generic feature of the potential structure suggests that instantons do not overlap 
in the spatial four dimensions although they look overlapped in the spatial three dimensions. 

Appendix G 

One-boson-exchange potential revisited 



In this appendix, we rederive the one-boson-exchange potential of §3 by summing up an 
infinite number of one-boson-exchange diagrams explicitly, in the standard field-theoretical 
computation. A key ingredient for this computation is the nucleon-nucleon-meson cubic 
couplings obtained in Ref. I8|),*3 It is found that the results are in total agreement with 
those derived in §3 See §4.3 of Ref. E]) for the definition of the couplings and the effective 
Lagrangian that we use to obtain the Feynman rule. 

G.l. Pion exchange 

The Feynman rule for the Yukawa coupling among a pion, nucleon N and nucleon N 
reads 

ig-KNNl^T"" (isotriplet sector) , idnNNl^T^ (isosinglet sector) . (G-1) 

Consider the scattering process where two initial nucleons with (pi,si,/i) and {p2,S2, I2) 
scatter to the final state composed of the two nucleons labeled as (p'j^, s'^^, ![) and (pg, Sj, I2) 
by exchanging a single pion. Here, Pi,P2,p'i, and are the on-shell momenta with the 
nucleon mass given by m^. si,S2,s'i, and specify the third components of the spin of 
the nucleons, and Ii,l2,I[, and stand for the third components of the isospin. It turns 
out that the scattering amplitudes due to the isotriplet and isosinglet pseudo-scalar meson 
exchange are given by 



X U{p[, s'i)75m(Pi, Si) — ^ U{P2, S2)J5U{P2, S2) , 



*^ In Ref. [H]), the cubic couplings involving excited baryons with / = J = 1/2 are also calculated. The 
extension of the computation of the nucleon-nucleon potential in this appendix to such excited states is 
straightforward . 
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Mfi)=Alf(^) ^ , (G-2) 

Sir JV]V-S>?7r AT JV, t"-!>tO 

respectively. Here, k = pi — p[ = P2 — P2 and 

rh = X^'''r^X^'\ (G-3) 

with x^^^^^^^ = (1)0)^ and x^^^'^^"^^ = (0)1)"^ being the isospin wavefunctions. The same 
notation will be used for the spin matrices 0"^/^,. For the definition of the Dirac spinors, see 
Appendix B.2 in Ref. [H]). Note also that we regard the pion as being massive with 7^ 
for the moment although the pion is massless in our model. 

In the large Nc and large A limit, the nucleon mass scales as 0{XNc) so that the 
nonrelativistic approximation is valid by considering the momenta to be of order one. Then 

Ei = E[ = E2 = E'^c^mB , 

~ -^{p? - p?) = 0{m~r}) , k^^P . (G-4) 

Furthermore, it can be shown that 

u{Pi, s;)75m(pi. Si) ~ {pi - p[)a . (G-5) 

Hence, 

>lf ^ H^mB? in ■ r,) (k . a,) (k ■ a,) , (G-G) 

[2mBY V /V ^ k^ + ml 

and a similar expression holds for the U{1) part. Here, in abbreviation, 

n = Ti'^i, , ^^ = ^s[s, ■ {i = l, 2) (G-7) 

To estimate the order of the amplitudes in A and N^, recall that the Yukawa couplings are 
given in Ref. 18|) as 

with f1 = (4/7r)/t. Here, ( ) denotes the expectation value with respect to the nucleon 
wavefunction given in §21 This implies that 

9nNN = 0{X'/' N'J') , g^MN = OiX-'/' iVi/2) , (G-9) 

showing that the SU(2) part dominates the U{1) part. 
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The effective action of the nucleons is defined to reproduce the amphtude computed 
above. In particular, the effective potential should be equated with 

1 



yirNN 



(G-10) 



mt 



(2mB)2 

Note that the overall factor (2772^)^ is removed from flG-6p because this comes from the 
wavefunctions assigned to the four external lines. By Fourier-transforming this, we obtain 

d^k 



Using the formula 



d^k 



ik-x 



(2vr)' 



1 



1 e 



(27r) 



)3 p + m2 47r r 
with r = and also for any function g{r), 



1 



1 



ai ■ V ■ V g{r) = - ((xl ■ a^) W'g{r) + -5i2 d^g - -Org 



with 



S 



(cTi • r ) (0-2 ■ r) 



12 



(G-11) 



(0-12) 



(G-13) 



(G.14) 



being the tensor operator, we find 



VJx) 



1 



ri ■ r2 



S 



12" 



47r (2mB)2 
Using the Goldberger-Treiman relation 



ml 



■ + — + ^ I + 



ml 



(o-i ■ 0-2) 



9a 



f-n-gTV 



NN 



TUB 



(G-15) 



(G-16) 



we find that the one-pion-exchange potential flG-lSp agrees with the expression (15-281) used 
in §3 As discussed in §Sl the central force vanishes when m.,^ = and the potential flG-lSp 
reproduces (l5-26p . and the n = component of the tensor force fl5-37p . For 7^ 0, it is 
standard in the literature to define the coupling 



4:71 \2mB J 



with which 



V^{x) = y(fi • T2) 



+ (o-i ■ 0-2 



(G-l?) 



(G-18) 
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G.2. Axial-vector meson exchange 

The Feynman rule for the nucleon-nucleon-axial-vector-meson cubic couphngs is 

da^NN "1151^ — (isotriplet sector) , da^NN ^151^ - (isosinglet sector) . (G-19) 

Here, a"' {n = 1, 2, ■ ■ ■ ) is the axial- vector meson associated with the wavefunction tp2n with 
the mass squared given by Asn- The propagator for a massive (axial-)vector boson of mass 
m is given by 

From these, the amplitude of the two nucleons exchanging an isotriplet axial-vector meson, 
summed over the species of the exchanged mesons, becomes 



4 in ■ nj 



E fc^ff^ U'^ + T^) Wi' ^i)757^^(Pi, ^i)) W2, 4)757'^n(p2, ^2)) , (G-21) 

n>l "T 2n V 271 / 

and we obtain a similar expression for the isosinglet case. Note that the cubic couphngs are 
computed in Ref. [8]) as 

dar^NN = \ {dzij2n{Z)) , Qa^NN = {p'^) {dz1p2n{Z)) . (G-22) 

This shows that 

and therefore the isosinglet sector is negligible compared with the isotriplet sector, as in the 
pion exchange case. 

In the nonrelativistic limit, where 

u{p', s'h^^Mp, s) = 0{m-^') , u{p', s'h,i'u{p, s) = %g{,^ + 0{mf) , (G-24) 

the scattering amplitude is dominated by the spatial component of the axial- vector fields so 
that 

Alf ^{2msf \ (n ■ f,) ±^ [(a, ■ + Y- (^^ ■ ^0 ■ 

ri>l ~^ 2n L 2n 

= -(2mB)'vif^(') . (G-25) 
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This yields the effective potential due to the axial-vector-meson exchange: 



X) 



1 



n>l 



(2^ 



Ak-x 



lai ■ (72 



k ■ Gi] [k ■ a2 



k^ + \2n ^2„ fc2 + A 



n>l 



NN 



s 



12 



1 + V-^2n'^ + ^A2„r^ 



1 

3' 



2n 



3 If^i ■ (^2j 



where 



^c"^ = ^ (^1 ■ ^2) (<?! ■ ^2)^glr^NNy2n{r) , 

n>l 



+ A2n >^2n('^ 



(G-26) 

(G-27) 
(G-28) 



Using (IG-22p . it is easy to show that this agrees with the n = 2, 4, 6, ■ ■ ■ components of fl5-36p 
and ([53ZD- 

G.3. Vector meson exchange 

The Feynman rule states that for the nucleon-nucleon-f " cubic couplings, we assign 

1 



gv"NN 7 ~ 



2m 



B 



Qv^NN 7 ~ 



2m 



(G-29) 



B 



for the isotriplet and isosinglet cases, respectively. Here, f " (n = 1,2, ■■■) is the vector 
meson associated with the wavefunction ip2n-i, whose mass squared is equal to A2n-i, and k 
is the momentum of the vector meson flowing outwards from the vertex. It follows that the 
two-nucleon scattering amplitude due to the exchange of an isotriplet vector meson, summed 
over the infinite tower of the vector meson species, is given by 

-1 



A1^^(2) = v^V^V^V^ — ■ r2) 



n>l 



{p'l, s'l) ( gv^NN 1^ - ^^^^p ) "bi> -Si) 



k^ky 



A 



2n-l 



X U{P2, S2) ( gv^^NN 7 + 0- k^ ] U{p2, S2) 



(G-30) 



As before, a similar expression follows for the U{1) part. 
In the nonrelativistic limit, we have 



u{p', s')7° u{p, s) = -i5s's + 0{m^^) 
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u{p', s')7^ u{p, s) = {p + p')j 6ss' + itjiaip - P')i + C'(m/) , 

u{p', s')cr°^ u{p, s) = (p - p')j Sss' + iejia{p + p')i (y^s's + , 

ZTTIb l -I 

u{p', s')a^'' u{p, s) = -e^^'' + 0(m^') • (G-3i; 
Furthermore, we note that the cubic couphng constants obtained in Ref. [8]) are given by 



1 (^2n-l(Z)) 



(G-32) 



which imply 

g^.NN = 0{X-'/'Ny') , h^.^M = 0{X-'/'N'J') . (G-33) 

Then, among the spinor bihnear forms appearing in the amphtudes, the leading ones for 
large A and large Nc are 



h 



u{p, s') ( g^nNN 7^ - 'T^^ - P)p ] u{p, s) 



2m 



B 



h. 



v"NN 

2171 R 



e,ia{P~P)Ks + --- , (G-34) 



u 



{P, s') I dv^NN 7° - 2^^^ cr'^''{p - P)p I u{p, s) ~ -igv"NN Ss's H • (G-35) 



This shows that for the isotriplet vector mesons, the spatial components dominate the am- 
plitude, while for the isosinglet vector mesons, the time component does. Consequently, 
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Fourier-transforming the effective potentials gives 
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Again, this is in agreement with the n = 1, 3, 5, ■ 



components of fl5-36p and f l5-37p . 
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Note added : When preparing this paper, we became aware of Ref. [57|) . which overlaps 
partly with our strategy. 

Note added in the second version: The added appendix [G] has some overlaps with 
Ref. EHJ, which appeared when we were preparing the revised version. 



56 



